Confinement of light quarks in QCD at finite density. 
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Abstract 

Cy . 

The 4q effective Lagrangian and the gap equation are derived for light quarks 
in the confinement phase of QCD. The modification of the confining string due to 
finite quark density (chemical quark potential /i) is observed. As a surprising result 
in a multiquark system with a common string junction an attractive well appears of 
radius \ija and of an average depth equal to fi. Possible implications for the density 
^NJ ■ phase transition are discussed. 
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§ ! 1 Introduction 

p , 

Recently the high-density effects in QCD attracted a lot of attention because of the pos- 
sible density phase transition pQ, and interesting accompanying phenomena like Color 
Superconductivity (CS) [2j (for earlier papers see JH|). On the fundamental level the 
physical expectation of any phase transition may be connected to the possible reconstruc- 
tion of the vacuum and for that one needs that the energy density to be of the order 
of the vacuum energy density e cr ~ e vac ~ — yiV c ^:((F" ) 2 ), e cr ~ 1 GeV/fm 3 which 
provides the drastic change in the vacuum structure and may cause the phase transition. 
In case of temperature phase transition the corresponding energy density e is indeed of 
the order of e cr , and as it was first argued in jij and later measured on the lattice jS], 
the QCD vacuum is strongly transformed in a way, that most part of colorelectric fields 
evaporate above T c . In this way T c was calculated through e cr and found in good agree- 
ment with lattice data |3]. Also, it was shown that in QCD in the framework of hadron 
resonance gas and making use of the low-energy theorems at T ^ [S] the quark conden- 
sate and one half (colorelectric component) of gluon condensate evaporate at the same 
temperature [Zj, which corresponds to the temperature of quark-hadron phase transition. 
Besides application of the effective dilaton Lagrangian to gluodynamics jH] and to QCD 
with light and heavy hadrons jHj permitted to conclude that similar reconstruction of the 
nonperturbative gluon vacuum takes place at finite temperature. 

A similar arguments for the density phase transition would imply that at the baryon 
density of ~ 1 nucleon/fm 3 , i.e. 3 -r 6 times higher than the standard nuclear density, the 
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vacuum can be reconstructed in such a way, that part of fields, e.g. the colorelectric fields 
responsible for confinement, disappear above critical density. 

A completely different route was trodden by the groups who studied the phenomenon 
of color superconductivity ,2J. Here main argument was the model study of CS in the 
formalism of the NJL or instanton model (see ^Hj for a recent review) and the practically 
important question of coexistence of quark and nuclear matter was also studied lately in 
the unified NJL approach ^T] . In this strategy the question of the vacuum reconstruction 
is not addressed, however it is usually assumed that confinement disappears in the course 
of establishing of the CS dynamics. On the other hand the possible importance of vacuum 
fields for the density transition and CS was suggested in [T2], but no detailed theory was 
presented there. It is the purpose of the present paper to start the investigation of the 
role of density on the vacuum fields in general and confinement in particular. 

In this paper we ask ourselves a short and simple question: how the nonzero baryonic 
chemical potential /x acts on the confinement of light quarks, and come to the unexpected 
answer, that the confining string of light quarks is destroyed gradually by /x in such a way, 
that the one part of string, near the string junction is eaten by the nonzero /x, while at 
distant r the string survives. The plan of the paper is as follows. In section 2 the effective 
Lagrangian is derived from the QCD Lagrangian with nonzero /x and their solution is 
discussed in section 3. Section 4 is devoted to the physical implications of results and 
prospectives. 



2 Derivation of Effective Aq Lagrangian 

One starts with the QCD partition function in presence of quark chemical potential /x in 
the Euclidean space-time, and we begin with the zero temperature, T = 0. 

Z = J DADlpDlp + e~ So ^ + $ S '4> + {id+im-i^ i +gA) fi/jd 4 x ^ 

where Sq(A) = \ §(F^ y (x)) 2 d A x, m is the current quark mass (mass matrix m in SU(3)), 
and the quark operator *ip aa (x) has flavor index a(f = 1, ■■■rif), color index a(a = 1, ...N c ) 
and Lorenz bispinor index a(a = 1, 2, 3, 4), and we use the contour gauge [T3] to express 
A p (x) in terms of F pv . One has for the contour z^s, x) starting at point x and ending at 

Y = z(0,x) 

A,(x) = / 1 ^ ^f; X) ^; X) ^(^(^)) - f Y dT» vp {z)F vp {z). (2) 

Integrating out the gluonic fields A p (x), one obtains 

Z = J Di/iDi/j+el f ^ + ( i 9+i^-i^i) f i'd 4 x e L^l }L +L^ ) QL +... ^ 

where the EQL proportional to ((A 11 )} is denoted by L^q L , 

L% L = C / d'xd'y tfaix) t^pix) ^+(y) 9 My)(^\x)A^\y)h^J (4) 
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Average of gluonic fields can be computed using (J2J) as (see [13] for details of derivation) 

g 2 (A^(x)A^(y)) = ^ f dma^u) f dv k a u (v)D^ H \u - v){8^5 ik - 5 iv 5 klx ), (5) 

l\ c Ja Jo 

where D <yE,H \x) is the correlator (Ei(x)Ei(0)) or (Hi(x)Hi(0)). As it was argued in [T3j 
the dominant contribution at large distances from the static antiquark is given by the 
color-electric fields, therefore at the first stage we shall write down explicitly L^q L (el) for 
this case, i.e. taking /i = v = 4. As a result one has[L4~] 

L% L (el) = ^ / / d 'y f ^(x) f A,(x) 9 ^{y) 9 ^(yh^HP E (x, y) (6) 
where J E (x, y) is 

J E (x,y) = [ X dm [ V d Vi D E (u-v), i = 1,2,3. (7) 
Jo Jo 

One can form bilinears = 9 Vv an d project using Fierz procedure given isospin 

and Lorentz structures, — > \I/( n ' fc )(a;, y). Here we consider only ip + ip bosonization. 
With the help of the standard bosonization trick (here J = jfJ E ) 

e -*i* = J (det jy/2 Dxexp[ _ x j x + mj x + %x j^ (8) 

Z = J DipDip + Dx exp Lq M l (9) 
one obtains the effective Quark-Meson Lagrangian (QML) 

l qml = J " d 4 x J d A y { f ip+ a (x)[(id + im - i^) a p5{x - y) + iM ( J^\x, y)\ 9 ^{y)~ 

-^-X (n ' fc) (x, y) J E (x, y) X {n ' k) (y, *)} (10) 
and the effective quark-mass operator is 

M% 9 \x,y) =J2x (n ' k) (x,y)OiXJj(x,y). (11) 

n,fc 

Here the operator O is a set of all irreducible combinations of Dirac matrixes. 

The QML in Eq.lflOj) Lq ML contains functions ^( n,fe ) which are integrated out in ([§]). 

and the standard way is to find x^'^ from the stationary point of Lq\, il . Limiting oneself 
to the scalar and pseudoscalar fields and using the nonlinear parametrization one can write 
for the operator M in (fTU|) 

M(x,y) = M s (x,y)U{x,y),U = exp(i>y 5 (f)),(f)(x,y) = <p f (x,y)t f . (12) 

After integrating out the quark fields one obtains the ECL in the form 

Lf CL {M s ,i) = ~2n f N c {J E {x,y))- 1 Ml{x,y)+N c tr\og[{id + im-i^ lA )l + iM s lJ]. (13) 
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The stationary point equations 5 ^ L = = at <p — 4>o, M s = immediately 

show that 0o = and Mj°) satisfies nonlinear equation 

iMf (x,y) = 4trSJ E {x,y) = {^S l4 )J E {x,y), S{x,y) = -[id + im - + %M s U]-\. 

(14) 

This equation plays the role of the gap equation and is the main point of our further 
investigation. For /i = this was done in ^1] and in the next section we find how results 
of |T3] are modified by the nonzero \i. 



3 The confining string at nonzero /i 

Our basic equations ©,© are nonlocal in time because of the integral over dx 4 dy 4 in 
(10). This nonlocality and the parameter which it governs can be handled most easily, 
when one uses instead of M(z, z'), S(z, z') the Fourier transforms. 

S{z 4 - 4, z, z') = [ e^-^S(p„ z, z')^ (15) 
and the same for M(z, z'). Then from (|14|) one obtains a system of equations 

(p4 — id z — im + ifX'j 4 )S(p 4 , z, w) — i J M(p 4 , z, z')S(p 4 , z', w)dz' = 5^ 3 \z — w) (16) 

To simplify matter, one assumes for D E (x) the Gaussian form, D E (x) = 
-D(O) exp ( — JfyJ- Then for M(j> 4 , z, w) one has 



*M(p 4 ,z, w) = 2^T g [ d A e -^-v',?n 

J 2tc 



x 



x [J E (z,w) l4 S(p' 4 ,z,w) l4 ] (17) 

where J E is defined in © and we have factored out the time-dependent exponent, using 
the Gaussian representation of D{u). 

All dependence of M onp 4 as can be seen in (|T7|) is due to the factor exp[— (p 4 — p' A ) 2 Tg] 
and disappears in the limit when T g goes to zero, while the string tension o ~ D(0)T 2 
is kept fixed. This limit can be called the string limit of QCD, and we shall study its 
consequences for equations (fl6 ^1 .(fT7 jl in this section. 

So in the string limit, with M independent of p 4 , let us consider the Hermitian Hamil- 
tonian 

H^j n = {-j-— +(3m- /i)^n(z) + P J M(p 4 = 0, z, z')^ n {z')d z z' = e ft (/x)V>„(z) (18) 
with eigenf unctions ip n satisfying usual orthonormality condition 

4>+{x)ip m {x)d 3 x = 6 nm , 
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From (jTgj) it is clear, that one can redefine e n {fi) + fi = e n , and e n and ^ n do not depend 
on /i. Therefore in all subsequent formulas one can use the same equations as in ^3], but 
with the replacement e n — > e n — fi. In particular, the Green's function S can be expressed 

as 

5(p 4 ,x,y) = 2^ 7 r— (19) 

n P474 - H^n - /i)74 

Inserting (fTTJJ) into (|T7jl one has integrals of the type: 

/OO (|p' g — (P4~ P4) 2 Tg ^ 

-oo 2vr (p' 4 7 4 - t(e n - /i)7 4 ) 2 

Note, however, that the result depends on the boundary conditions. If, e.g., one imposes 
the causality-type boundary condition, then one obtains 

dp' A e< hi _ f i<y 4 e- £h *6(e n - fi), h 4 > 

2vr 7 4 (p' 4 - i(e n - ^)) ~ \ ~i^0(fi - e n )e e?t4 , /i 4 < 

We are thus led to the following expression for M in the string limit 

M(p 4 = 0, z, w) = y/^T g J E (z, w) 74 A(z, w) (21) 

where the definition is used 

A ( z > w ) = ^n{z)sign{e n - /i)^(w) (22) 

n 

Let us disregard for the moment the possible appearance in M of the vector component 
(proportional to 7^,/U = 1,2,3,4) and concentrate on the scalar contribution only, since 
that is responsible for CSB and confinement. Then one can look for solutions of the Dirac 
equation (|18jl in the following form 

Mrl = -(f/^ lM ) (23) 

where V = 2j — I, and introducing the parameter k(J, I) = (j + ^)sign(l—j), and replacing 
M by a local operator (the generalization to the nonlocal case is straightforward but 
cumbersome, for a possible change in the nonlocal case see [Hj). We obtain a system of 
equations 

J ^ + *G n -(e n -fi + m + M(r))F n = 
{ ^-*F n + (e n -fi-m- M(r))G n = 

Eq.(j23J) possesses a symmetry (e n — fi, G n , F n ,K) <-> (fi — e n ,F n ,G n ,—n) which means that 
for any solution of the form ()22j) corresponding to the eigenvalue e n — ft, there is another 
solution of the form 

= 1 f feMsT ) < 25 > 

corresponding to the eigenvalue (fi — e n ). 
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Therefore the difference, which enters (J22|) can computed in terms of F n , G n as follows 

A(z, w) = A (z,w) - AA(z,w) (26) 
where A is the value of A for \i = 0, i.e. the same as in ^3], while AA is defined as 

AA(z,w)=2 ]T V«(z)V+(w). (27) 

0<£ n </i 

Using decomposition (J23j) one can write AA as 

where we have denoted Q = fljiM, ^' = Qji'M, and we disregard nondiagonal part of A. 

At this point one can follow the relativistic WKB method for Dirac equation |T5] 
applied to calculation of A in [14J in case of \x = 0. The classically available region for 
ip n (r) with energy e n = e is r min < r < r max , where r maXjmin = s 2 ±Ve-4^^ ; an d the 
summation over n in (}2"Tj) transforms into integration over de, with the lower limit (for a 
given r) e m i n = or. In this way one has for the upper diagonal element in ()28|) . 

2o /-At/o-r r + 1 / 

AA(+, +) = -=-5(1 - cos^ rr / dr . = cos{aVr 2 - l)9U - or) (29) 
n 2 r Ji yr 2 — 1 

with a = ar\r — r'\, and we keep r m r' everywhere except for a, since for large a (when 
r is far from r') both A and AA fast decrease. 

In a similar way for the lower diagonal element in (|28|) one has 



AA(— , -) = —5(1 - cos r r') / dr - cos(aV^ T ^l)9(fi - or) (30) 

7t 2 r Ji yr 2 — 1 

and taking A in (|25|) from JT3j the resulting form for A (|26|) is 

/3<7 , . i" 00 drcos(a\/ 



a/ /\ i, per r/1 n \ f°° "Tcosla 

A(r, r ) = fJA scalar + lA„ ector = -=-d(l - cos 6> rr /) / = 

7rr j/i/ot vr 



r 2 



2cr „ , /"^/' jr rdr 



- 1—5(1 -cos Orr) cos(ay / r 2 - 1). (31) 

7r z r ji yr 2 — 1 

Here (3 = 74, and 1 is the unit Dirac matrix, which means, that the second term on 
the r.h.s. of (|3*T|) contributes to the vector part of the resulting mass operator (J3*|). while 
the first term contributes to the scalar part. One should take into account, that in the 
first term, 

°" r/, n s f°° dr cos(av / r 2 - 1) , . 

Asca/ar = ^1 " COS ^ / )^== 1 32 

vr 2 r A min ( M ) Vr 2 - 1 

T min(n) = /i/crr for fi > or and 1 otherwise, so that for large r, r 3> M , one has the 
standard //-independent value 

Kcaiar(r ~r' > fi/a) = —K (a)5(l - cos0 rr ') (33) 

7r^r 



where we have used relations for the McDonald function K 

= / /i-T-2 ' / daK °(a) = o- (34) 
•'O v 1 + x 2 Jo 2 

One can check that at large r, r'(r ~ r' > /x/cr) A scaiar « A^°|. is a smeared 5 
-function, 

AS>(r,r')dV = l. (35) 



However for large /i, /i ^> y^o 7 , A sca i ar is different from A sca i ar (fi = 0), and for r ~ r' < fi/a 
one has approximately 

Kcalar(r, r') = - /(r, r')0(// - <«0 (36) 

where 

f{ r , r ) = / s 5 =/ , , 9 cosaA, (37) 



and A = ^) 2 -l. 

Let us now turn to the vector part of interactions, A„ ector , 

A„ ector . = ^-^(l - cos6> rr ') / - cos(aVr 2 - 1)6 (n - err) 

7T z r Ji Vr 2 — 1 



On- sm(m/f JL> ) -1) 
1^(1 - cosg rr VU; -. (38) 

Returning back to Eq. (|2ip one can deduce, that 

M(p 4 = 0, r, r) = \fnT g J E {r, r') [A SCQi (r, r') + 7 4 A 1;ector ] = M smi + 74M„ ect . (39) 
Taking into account, that at r, r' ^> T 5 and for the Gaussian D E (x) one has from (J7| 

J B (r, r') = ^— /2T 9V ^F J D(0) min(r, r') (40) 

one has for M at r, r' ^> T g and for r = r' 

M scal (r,r') = ar(8®(r,r')-Z(r,r')) (41) 
where <5( 3 )(r, r') = A^ (r, r') and 

£(r, r') = 5 M = ^-(l - cos0rr')/(r, r')0(/i - err) (42) 



7r^r 



and f(r,r') is defined in (J37J). 
For M vect one has, using 



M vect (r, r') = -2ar^ M (r, r')0(/i - <rr) (43) 



where 

sm(a\/f^) 2 - 1) 
^(r,r') = —5(1 - cos6 rr ) . (44) 

Note, that one should actually symmetrize all these expressions, e.g. - — > -7= etc., 
but we always are in the regime, where r m r' . To estimate the magnitude of nonlocal 
kernel M sca i(r, r') and M vect (r, r') it is convenient to introduce in (|18jl the local limit of 
the mass operator, namely 

A»,,« d (r)=/dVA^,,„„ i (r,r'). (45) 

Exploiting the equalities 

J dW 3) (r,r') = J dV£(r,r') = J d 3 rV /t (r, r') = 1 (46) 

one arrives at the expressions 

M scal (r) = ar9{ar - fx), M vect = -2ar9(fi - err). (47) 

In the next section we shall discuss approximations made in deriving (JUJ, (|4*3|) . (jUJ) 
and physical implications of these results. 



4 Discussion of results 

Results of the previous section Eqs. (|41|) . ()4llj1 . (|47jl . can be formulated as follows. The 
relativistic WKB analysis leads to the /^-dependent modification of the confining string, 
where the piece [0,/x/cr] of the string near the origin of the string (situated at the heavy 
quark position in case of heavy-light quark, or at the string junction position in the case 
of baryons), is dissolved, and the linear confinement starts beyond the critical radius 
r C r = /V* 7 - Moreover, an attractive vector interaction appears in the same interval with 
the average magnitude (M vect ) ~ fi. 

These conclusions should be taken as qualitative. First of all, the WKB method is 
not a good approximation at small distances, and we have omitted exponentially damped 
part of ip n (z) in the spectrum, therefore the inner part of the string is to some extent 
delocalized (see [Hj for details) and smoothed. 

Secondly, we have not taken into account a possible modification and destruction of 
the vacuum due to the influence of high density quark matter, which might decrease o or 
cancel the string completely (as it is happens in the thermal phase transition jlj). The 
phenomenon of this kind was observed on the lattice jTH] where deconfinement tempera- 
ture decreased under the influence of applied external Abelian field. 

If however, no density induced vacuum deconstruction takes place, then the resulting 
physical picture according to Eqs. (|47jl . is the net decreasing of confinement in the inner 
region of some ensemble of quarks, and appearance of attractive vector potential of the 
order of fi acting on each quark. This may cause creation of deconfined bubbles consisting 
of 3n quarks, n = 2, 3, .. in the midst of the nuclear medium, and dynamically is similar to 
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the 3nq bag formation, which was studied before in the framework of the Quark Compound 
Bag model ^7j. Note, however, that bag boundary conditions might be strongly modified 
as compared to the standard MIT bag model. A quantitative analysis of this situation 
needs a more accurate analysis of the 3n quark system using nonlinear equations for the 
3nq Green's function, generalizing Eq. 1)14)1 . 

The formation of these high-density 3nq bubbles may be connected with the expla- 
nation so-called cumulative effects in the hadron-nucleus(and nucleus-nucleus) collisions, 
for an example of this discussion see ^H] and refs. therein. 
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